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Stopping  powers  derivable  from  Lindhard’s  dielectric  theory  of  stopping  for  a  charged  .*■ 

particle  in  a  uniform  electron  gas  are  calculated  as  a  continuous  function  of  electron  U.**' 

density  for  a  wide  variety  of  projectile  energies  ranging  from  102  to  104  keV/amu.  f  "  I 

Results  are  presented  for  electron  densities  varying  from  1022  to  1032  electrons/cm3,  a  ■S'® 
range  which  includes  electron  densities  generally  found  for  the  energy  bands  of  solids  APR 

and  for  the  atomic  levels  of  atoms.  <  — ,  W t 
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INTRODUCTION 

Lindhard'  presented  an  elegant  dielectric  formulation 
of  stopping  theory  for  a  charged  particle  penetrating  a  free- 
electron  gas.  His  work  probably  stands  as  the  most  complete 
“first  principles”  description  of  energy-loss  theory  to  date. 

The  original  work  of  Lindhard  was  later  expanded  by 
Lindhard  and  Wintheri  (L-W)  who  obtained  analytic  forms 
for  the  stopping  power  in  the  limit  of  high  and  low  projectile 
velocity  relative  to  the  electron  gas  Fermi  velocity.  This  ex¬ 
pansion  into  analytic  expressions  was  necessary  to  avoid  the 
significant  computational  difficulties  implicit  in  the  original 
Lindhard  formalism.  They  also  presented  results  obtained 
from  a  direct  numerical  integration  of  the  Lindhard  stop¬ 
ping  equation  for  six  different  electron  densities  spanning 
several  decades  in  value. 

In  this  paper,  results  of  the  numerical  integration  of 
Lindhard's  original  stopping  equations  are  presented  in  de¬ 
tail.  Further,  polynomial  fits  of  the  results  are  presented  so 
that  the  Lindhard  stopping  power  can  be  obtained  directly 
for  a  given  electron  charge  density.  These  results  are  present¬ 
ed  as  a  function  of  electron  density  so  as  to  facilitate  the 
calculation  of  stopping  powers  with  the  local  density  ap¬ 
proximation  (an  approximation  in  which  each  volume  ele¬ 
ment  of  a  target  is  considered  to  be  an  electron  gas  of  uniform 
density;  then  the  total  electronic  stopping  power  is  obtained 
by  an  appropriate  averaging  over  the  spatial  variation  in  tar¬ 
get  electron  charge  distribution). 

LINHARD'S  THEORY  OF  STOPPING  POWER 

For  an  ion  of  charge  ze  moving  with  velocity  v  in  a 
medium  of  uniform  density  p,  the  energy  loss  due  to  electron 
excitation  can  conveniently  be  written  in  the  form 


electron.  In  the  dielectric  formalism,  L  is  written  as 

L  =  —I—  r^-  f  a,  dco  [e  '(k,a>)  -  1  ],  (2) 

mogJo  k  J ~kv 

where  o>0  is  the  plasma  frequency 

...2  _  n.\ 


and  e(k,co)  is  the  wave-number-  and  frequency-dependent 
longitudinal  dielectric  constant. 

Lindhard  obtained  the  dielectric  constant  for  a  free- 
electron  gas  within  first-order  perturbation  theory  to  be 

e(k,w) 

,  .  8 irme2^„,„J  1 


-I^(£n)( 


dE  4W  re2  V  ..  , 

-  — =  -H —  pL  (p,v), 
dx  m  \  v  / 


where  L  is  called  the  stopping  number  and  m  is  the  mass  of 


k 2  +  2k-k„  -  (2m/fP)(a)  +  iy) 


k 2  -  2k-k„  +  (2 m/W)(p)  +  iy) ) 

where  En  and  k„  represent  the  energy  and  wave  vector  for 
an  electron  in  the  nth  state  and  y  is  a  small  positive  value.  For 
a  free-electron  gas  at  zero  temperature, 

F(En)=\  for  En<EF 

=  0  for  E„  >  Ef  (5) 

and  the  summation  in  Eq  (4)  can  be  evaluated  to  get 

e(uj)  =  1  +  (xVz2)  L /!(«>*)  +  ] ,  (6) 


f  (M  =  -7  +  7"[  1  -  (*  -  «021  In  - — 

2  8z  z  —  u—l 


+  f[l-(r  +  «)2]  ln*  +  M  +  |-  ,  (7a) 

8r  z  +  u—  1 
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TABLE  I.  Tabulation  of  the  A,  coefficients  in  Eq.  (16)  for  Kf,  where  V,  is  the  Fermi  velocity  defined  by  Eq.  (9),  and  for  various  projectile  energies. 


i 


[  Energy 

I  (MeV) 

AaX  10  1 

/I.XlO" 

/I,Xl0 

/LX10* 

A.x  101 

I  0.1 

-  5.229 

4.526 

-  1.078 

9.599 

-  3.268 

||  0.2 

-  8.132 

6.293 

-  1.467 

13.39 

-  4.651 

i  0.3 

-  11.96 

8.574 

-  1.968 

18  26 

-  6.425 

I  0.4 

-  10.36 

7.641 

-  1.759 

16.19 

-5.66) 

0.5 

-  29.12 

18.67 

-4.178 

39.73 

-  14.23 

0.6 

-  16.1 

11.01 

-  2.491 

23.26 

-  8.215 

0.7 

-  24.54 

15.94 

-  3.567 

33.66 

-  11.97 

!  0.8 

-  46.47 

28.77 

-  6.370 

60.82 

-  21.82 

1.0 

-  23.79 

15.48 

-  3.453 

32.46 

-  11.50 

2.0 

-93.14 

55.58 

-  12.12 

115.5 

-41.31 

3.0 

-  107.1 

63.42 

-  13.76 

130.7 

-  46.56 

4.0 

-  76.02 

45.32 

-9.815 

92.63 

-  32.82 

5.0 

-  168.1 

98.06 

-21.12 

200.2 

-  71.13 

6.0 

-83.91 

49.65 

-  10.70 

100.6 

-  35.51 

7.0 

-  154.7 

90.02 

-  19.31 

182.2 

-  64.43 

8.0 

-  229.9 

132.9 

-  28.45 

268.7 

-9510 

10.0 

-  144.3 

83.73 

-  17.89 

168.0 

-  59.15 

and 

i iru ,  for  z  +  u  <  1 

=  (jr/8z)[l -(z  —  u)2]  for  \z  -  u\ <  1  <z  +  u 

=  0  for  |z  —  u  \  5  1.  (7b) 

Here  the  variables  k  and  to  have  been  replaced  by  reduced 
variables 


X1  =  v0/m)„ 


(10) 


where  v0  =  ei/fi,  the  Bohr  velocity.  When  the  dielectric  con¬ 
stant  of  Eq.  (6)  is  substituted  into  Eq.  (2),  one  obtains 


IT  Jo  Jo 

v _ 2f(uj) 


[*2  +  x  2f(u^)V  +  [x'Muj)  Y 


(11) 


z  =  k  /2k f  and  u  =  <o/kvf , 

where  kF  and  vF  are  defined  by 

1  ,  ffkl  fp 

Ef  =  -f mv'F=— i  = 

z  zm  zm 

In  addition,  the  parameter  in  Eq.  (6)  is  defined  by 


(8) 


(9) 


This  expression  for  L  in  Eq.  (1 1)  is  the  Lindhard  stop¬ 
ping  number.  For  simplicity  in  integrating  the  z  dependence 
in  Eq.  (1 1),  we  let  z  =  u  +  a  for  a  fixed  u  and  factor  out 
(l/}'2)  to  get 

L  =  f  u  du(  (u  +  a)yF(u,a)  da,  (12) 

TTX‘Jo  J  u 


TABLE  II.  Tabulation  of  the  A,  coefficients  in  Eq.  (16)  for  Vu,„  >  V,  ,  where  Vy  is  the  Fermi  velocity  defined  by  Eq.  (9),  and  for  various  projectile  energies. 


Energy 

(MeV) 

AoX  lO’ 

/LXlO‘ 

4jX  10"1 

A,X  10-4 

A.XlO-' 

/I,X10» 

A.X101 

^,X10* 

4,X  10’ 

0.1 

-244.8 

371.9 

-  247.0 

93.73 

-  2221 

336.7 

-  318.9 

172.5 

-407.9 

0.2 

-  247.3 

366.2 

-  237,2 

87.77 

-  2029 

299.9 

-  277.1 

146.2 

0.3 

-  158.1 

232.5 

-  149.5 

54.91 

-  1260 

184.8 

-  169.4 

88.71 

-203.1 

0.4 

-90.85 

132.9 

-  84.96 

31.03 

-  707.7 

103  3 

-94.10 

48.97 

-  111.4 

0.5 

-  85.82 

125.1 

-  79.69 

29.00 

-  659.1 

95.81 

-  87.00 

45.11 

-  102.3 

0.6 

—  41.18 

59.71 

-  37.84 

13.70 

-  309.5 

44.74 

-40.38 

20.81 

-46.89 

0.7 

-  52.28 

75.71 

-  47.93 

17.33 

-  391.3 

56.51 

-  50.96 

26.25 

-  59  09 

0.8 

-  33.78 

48.81 

-  30.83 

11.12 

-  250.3 

36.05 

-  32.42 

16.65 

-  37.37 

1.0 

0.8278 

-  1.550 

1.209 

-0.5216 

13.73 

-  2.272 

2.316 

-  1.333 

3.327 

2.0 

-  0.5526 

0.5599 

-  0.2033 

0.0192 

0.7829 

-  28.74 

4.151 

-  0.2934 

0.8381 

3.0 

2.114 

-  3.161 

2.063 

-  0.7676 

17.82 

-  2.643 

2.445 

-  1.290 

2.973 

4.0 

0.6460 

-  1.025 

0.7039 

-  0.2740 

6.618 

-  1.016 

0.9700 

-  0.5262 

1.243 

5.0 

5.541 

-  7  .965 

5.004 

-  1.795 

40.21 

-  5.758 

5.150 

-  2.629 

5.867 

6.0 

8,225 

-  11.76 

7.345 

-  2.620 

58.35 

-  8.310 

7.389 

-  3.751 

8.324 

7.0 

3.808 

-  5.453 

3.413 

-  1.219 

27.20 

-  3.880 

3.456 

-  1.757 

3  904 

8.0 

5.177 

-  7.381 

4.599 

-  1.636 

36.34 

-  5.160 

4.575 

-  2.315 

5.121 

10 

6.605 

-9.366 

5.804 

-  2.053 

45.35 

-6.404 

5.646 

-  2.842 

6250 

5080 
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LINDHARD  STOPPING 


FIG.  I .  Plot  of  variation  in  stopping  number  {Eq.  ( 1 2)]  with  electron  density 

for  select  projectile  energies.  [(— )  1 00  keV/amu;  ( - )  1000keV/amu;(-  ■ 

--)  10  000  keV/amu). 


where 


F(u,a)  = 


Y!fi(u,a) 

g*(u,a)  +  \x2fAu,a)) 2 


and 


(13a) 


g(u,a)  =  (u  +  a)2  +  \2fAu,a).  (13b) 

In  Eq.  (1 3a), /,(u,fl)  is  defined  by  Eq.  (7a)  with  z  —  u  +  a  and 
Mu, a)  is  given  by 

Mu, a)  =  ±iru  for  2u  +  a  <  1  (14a) 


— — —  for  |a|a<l<2u+a  (14b) 

8  (u  +  a) 


=n8[a  —  a0(u)] 


da 


-  i 


Uu (M) 


05) 


where  a0(u)  is  the  solution  to  g(u, a)  =  0  for  a  fixed  u.  Thus, 
when  |a|  >  1,  the  use  of  Eq.  (15)  in  Eq.  (12)  reduces  the  dou¬ 
ble  integral  to  a  line  integral. 

We  present  the  result  of  the  numerical  integration  of 
Eq.  (12)  for  100  electron  density  values  in  the  interval  10”- 
10”  electrons/cm1  and  for  projectile  energies  in  the  range 
between  iOOand  10*  keV/amu.  The  results  are  presented  in  a 
“least-squares”  polynomial  fit  of  stopping  number  with  elec¬ 
tron  density  in  the  form 

where  the  A,  and  n  are  given  in  Tables  I  and  II.  The  polyno¬ 
mial  fit  is  presented  in  two  parts,  one  for  vum  <,vF,  and  one  for 
yion  >vf,  where  vF  is  the  Fermi  velocity  defined  by  Eq.  (9). 
Typically  illustrations  of  the  variation  in  stopping  number 
with  electron  density  and  energy  can  be  found  in  Fig.  1. 

SUMMARY 

We  have  presented  the  results  of  the  numerical  integra¬ 
tion  of  Lindhard's  stopping  number  over  a  wide  range  of 
electron  densities  and  projectile  energies.  The  results  are  giv¬ 
en  in  the  form  of  a  polynomial  fit  of  stopping  number  with 
electron  density.  The  results  of  this  paper  will  subsequently 
be  used  to  apply  the  Lindhard  theory  of  stopping  to  “real” 
solids  via  the  use  of  the  local  density  approximation. 
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